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Abstract 

We present a number of results about (finite) multiple harmonic sums modulo 
a prime, which provide interesting parallels to known results about multiple zeta 
values (i.e., infinite multiple harmonic series). In particular, we prove a "duality" 
result for mod p harmonic sums similar to (but distinct from) that for multiple zeta 
values. We also exploit the Hopf algebra structure of the quasi-symmetric functions 
to do calculations with multiple harmonic sums mod p, and obtain, for each weight 
n < 9, a set of generators for the space of weight-n multiple harmonic sums mod p. 

1 Introduction 

In recent years there has been considerable interest in the multiple zeta values 

1 

n.i>n,2>--->nj.>l i ^ K 

which converge for ii > 1; see, e.g. [21 [131 [H [13 [HI [271 [21] • The multiple zeta values 
are limits of the finite harmonic sums 

A(^i,,i,,...,i,){n) = J2 -^TIT ^' (2) 

n>ni>n2>->nfc>l i ^ ft 

as n ^ oo. In fact, these finite sums themselves are of interest in physics [H [26]. Both 
types of quantities are indexed by the the exponent string (ii, . . . ,1^): we refer to k as 
the length and zi + ■ ■ ■ + as the weight. 

The multiple zeta values (1) form an algebra under multiplication. In fact, there are 
two distinct multiplications on the set of multiple zeta values, the harmonic (or "stufHe" ) 
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product and the shuffle product, coming from the representation of multiple zeta values 
as iterated series and as iterated integrals respectively (see §4 of [I?] for a discussion). 
Only the first of these products applies to the finite sums (2). 

In this paper we develop a theory of the mod p values {p a prime), of the finite harmonic 
sums (2) when n = p — 1. Our main tools are the structure of the algebra QSym of quasi- 
symmetric functions (which formalizes the harmonic product) and a mod p duality result 
(Theorem 4.7 below) which corresponds to the duality theorem for multiple zeta values 
(see, e.g.. Corollary 6.2 of [13]). The divisibility properties of finite harmonic sums are 
considered by Zhao [29j, and there is some overlap between his paper and this one. But 
in this paper we focus only on congruences mod p, while [29j also considers congruences 
mod higher powers of p. A brief account of the results of §4 of this paper appeared earlier 
in the last section of [16J. 

The theory of mod p finite harmonic sums has several interesting differences and 
similarities to the theory of multiple zeta values. While the values of (1) with length 
k = 1 are the classical zeta values C(^i)) it is well known that the mod p value of 

P-i ^ 

%)(p-l) = E-n 

is zero when p > ii + 1. But for A; = 2 we have the relation 

^{ti,i)iP - 1) = 5p-n-i mod p, 

with Bernoulli numbers (see Theorem 6.1 below). In the case of multiple zeta values, it 
appears that that multiple zeta values not expressible in terms of the classical zeta values 
occur first in weight 8 (e.g., ({6,2)); for the finite harmonic sums mod p, the first ones 
not expressible in terms of Bernoulli numbers appear to occur in weight 7 (see §7 below). 

As we have already mentioned, there is a duality result in both theories. If we code the 
exponent strings {ii, 22, • • • , ik) occurring in (1) and (2) by the monomial x^^~^y ■ ■ ■ x^'^~^y 
in noncommuting variables x and y, then the duality map for multiple zeta values is the 
anti-automorphism of Q,{x,y) that exchanges x and y, while that for the mod p finite 
harmonic sums is the automorphism ip of Z(x, y) with iIj{x) = x + y and ip{y) = —y. 

The dimension dn of the rational vector space of multiple zeta values of weight n was 
conjectured by Zagier [28] to be the Perrin numbers, i.e. di = 0, d2 = ds = 1, and 
dn = dn-2 + dn-3- The corresponding numbers for the finite harmonic sums (2) are the 
minimal numbers c„ of weight-n harmonic sums needed to generate all weight-n harmonic 
sums mod p for p > n + 1. From the calculations of §7 below the first few values of Cn 
seem to be the following. 
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(The value Cg = 2 is based on Theorem 7.5 below together with a relation conjectured 
by Zhao; see the remarks following Theorem 7.5). We do not have a conjecture about 
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the Cn for general n, but we do conjecture that all multiple harmonic sums Aj{p— 1) can 
be written mod p as sums of products of sums of the form — 1) (what we call 

"height one sums" ; see §5 below) with h even. 

This paper is organized as follows. In §2 we give an exposition of the algebra QSym 
of quasi-symmetric functions. In §3 we describe the Hopf algebra structure of QSym, 
particularly the antipode, and introduce some integral bases of QSym important in the 
sequel. In §4 we define bases Ai and Sj for the set of finite harmonic sums, and prove 
some basic results about the mod p case. In §5 we discuss a particular class of multiple 
harmonic sums, the height one sums. In §6 we obtain some results about Si{p — 1) mod 
p when / has short length. Finally, in §7 we apply the preceding results to find sets of 
products of height one sums that generate all multiple harmonic sums mod p through 
weight 9. 

This paper was written during a stay at the Max-Planck-Institut fiir Mathematik in 
Bonn, and the author thanks the directors of the Institut for their hospitality. The author 
also thanks Joe Buhler for his prompt response to a query about irregular primes. 

2 Harmonic Algebra and Quasi-Symmetric Functions 

Let S) be the underlying graded abelian group of the noncommutative polynomial algebra 

Z{x,y), where x and y both have degree 1. For a word w of 53, we refer to the total 
degree of w as its weight (denoted \w\) and its y-degree as its length (denoted i{w)). We 
can define a commutative multiplication * on by requiring that it distribute over the 
addition and that it satisfy the following axioms: 

HI. For any word w, l*w = w*l = w; 

H2. For any word w and integer n > 1, 

X * w = w * X = vox ] 

H3. For any words wi, W2 and integers p, g > 0, 

x^ywi * x'^yw2 = x^y{wi * x'^yw2) + x'^y{x^ywi * W2) + x^^'^^^y{w\ * W2). 

Note that axiom (H3) allows the *-product of any pair of words to be computed recursively, 
since each *-product on the right has fewer factors of y than the *-product on the left-hand 
side. Induction on length establishes the following. 

Theorem 2.1. The *-product is commutative and associative. 

We refer to S) together with its commutative multiplication * as the (integral) harmonic 
algebra (S), *). Let S^^ be the additive subgroup Z1 + Sjy of ij; it is evidently a subalgebra 
of (i^, *). Note that any word w E can be written in terms of the elements Zi — x^~^y, 
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and that the length i{w) is the number of factors Zi of w when expressed this way. We 
can rewrite the inductive rule (H3) for the ^-product as 

ZpWi * ZqW2 = Zp{wi * ZqW2) + Zq{z.pWi * W2) + Zp+q{Wi * W2). (3) 

Now for each positive integer n, define a homomorphism 0„ : ^ Z[ti, . . . ,tn] of 
graded abelian groups (where = 1 for all i) as follows. Let 0„(1) = 1 and 



■ ■ Z, 



E 



l<ri\<n2<--<nk<n 



for words of length k < n, and let (f){w) = for words of length greater than n; extend 
(pn additively to S^^. Because the rule ([3]) corresponds to multiplication of series, 0„ is an 
algebra homomorphism of [S)^, *) into Z[ti, . . . and (pn is evidently injective through 
degree n. For each m > n, there is a restriction map 



Pm,n '■ . . . , tm] Z[ti, . . . , t„j 



defined by 



The inverse limit 



Pm,n (ti 



ti, i < n 
0, i > n. 



<P = proj limZ[ti, . . . ,t„j 



is the subalgebra of Z[[ti, ^2, • • • ]] consisting of those formal power series of bounded degree. 
The (pn commute with the restriction maps, so they define a homomorphism : — 
Inside ^ is the algebra of symmetric functions 

Sym = proj lim Z [t i , . . . , t„] ^" 

n 

and also the algebra QSym of quasi-symmetric functions (first described in [S]). We say 
a formal series p G ^ is in QSym if the coefficient of t^^ ■ ■ ■ tf* in p is the same as the 
coefficient of t^l ' ' ''^^1 P whenever ii < i2 < ■ ■ ■ < ik and ji < j2 < ■ ■ ■ < jk- Evidently 
Sym C QSym. An integral basis for QSym is given by the monomial quasi-symmetric 
functions 

i\<i2<--<ik 

which are indexed by compositions (pi, . . . Since (p){zi^ ■ ■ ■ ZiJ = M(j^^ j^), we have 
the following result. 

Theorem 2.2. (p is an isomorphism of{S)^,*) onto QSym. 

As is well known, the algebra Sym of symmetric functions is generated by the ele- 
mentary symmetric functions Cj, as well as by the complete symmetric functions hi. The 



4 



power-sum symmetric functions Pi generate Sym (8>Q, but only generate a subalgebra of 
Sym over Z. It is easy to see that 4>{z\) = and 4>{zi) = Pi. 

A monomial quasi-symmetric function M/ is in Sym exactly when all parts of / are 
the same. Given / = {ii, . . . the symmetric group on k letters acts on I by 
a ■ I = (v(i), . . . , and the symmetrization 

of Mj is evidently in Sym. Hence the element (jlj) can be written as a sum of rational 
multiples of the power-sums Pi = M{iy But in fact (jlj) is a sum of products of the M(j) 
with integral coefficients that can be given explicitly as follows. 

Theorem 2.3. Let I = (ii,«2, • • • be a composition. Then 

E = E (-l)'=^'c(S)M(,oM(,,) ■ ■ ■ M(,,), 

o-eSfe partitions 2 = {Bi,...,Bi} of{l,...,k} 

where c(!B) = (card-Bi — 1)! • ■ ■ (card 5; — 1)! and bs = Xljes ^i- 

Proof. While the argument used to prove Theorem 2.2 of [T3j can be adapted to prove 
this result, we shall use the Mobius inversion formula as follows. The set Uk of partitions 
of the set {1, . . . , k}, with the partial order given by refinement, is a finite semimodular 
lattice. We also have, for any C = {Ci, . . . , Cp} G 11^, 

M(,,)M(,,) ■ ■ ■ M(,^) = E E ^-■ib„...Mh 

where bg = X]je-Bs '^^ ~ ^iec Then the Mobius inversion formula gives 

E ^-(^1--^.) = E /i(S,e)M(,,)---M(,,). (5) 

o-eSp 'B={Bi,...,Bi}^e 

The conclusion follows by taking C = {{1}, {2}, . . . , {k}} in equation ([5]) and noting that 
in this case e) = (-l)*^-'c(S) (see [21], Example 3.10.4). □ 

The rational algebra of quasi-symmetric functions QSym (8)Q was shown to be a poly- 
nomial algebra by Malvenuto and Reutenauer p2j. To describe the generators, we put 
an order on the words of by setting x < y and extending it lexicographically. A word 
is called Lyndon if it is smaller than any of its proper right factors, i.e. for w Lyndon we 
have w < V whenever w = uv for u,v 1. Let £ be the set of Lyndon words in Sj^. Then 
the following result was proved in [22j. 

Theorem 2.4. QSym(S)Q is a polynomial algebra on the generators (piw), where w (z L. 

The integral structure is more subtle. The result is again that the algebra is polyno- 
mial, but one has to use a different set of generators. Let L™-"'^ be the set of modified 
Lyndon words in i.e., whenever w = zf^ ■ ■ ■ zf^ G £j has a common factor d dividing 

all the exponents Pi, replace it with v'^, where v = zf^^^ ■ ■ ■ zf'"'^'^. The following result was 
first stated by Ditters, though the first correct proof seems to be due to Hazewinkel [12]. 

Theorem 2.5. QSym is a polynomial algebra on the generators (t){w), where w G , 
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3 Quasi- Symmetric Functions as a Hopf Algebra 



For definitions and basic results on Hopf algebras we refer the reader to [23], [25], and 
[20] . The algebra (i^^, *) = QSym has a Hopf algebra structure with coproduct A defined 
by 

n 

i=o 

and counit e with e{u) = for all elements u of positive degree. This extends the well- 
known Hopf algebra structure on the algebra Sym (as described in [?]), in which the 
elementary symmetric functions y*) and complete symmetric functions hi are di- 

vided powers, while the power sums Pi Zi) are primitive. The Hopf algebra (^^, *, A) 
is commutative but not co commutative. Its (graded) dual is the Hopf algebra of noncom- 
mutative symmetric functions as defined in [S]. 

Now QSym has various integral bases besides the Mj, also indexed by compositions. 
For compositions /, J, we say I refines J (denoted I >~ J) if J can be obtained from / by 
combining some of its parts. Then the fundamental quasi-symmetric functions are given 
by 

Fr = J2Mj, (7) 

and the essential quasi-symmetric functions are given by 

Ej = J2Mj. (8) 

It will be useful to have some additional notations for compositions. We adapt the 
notation used the previous section for words, so for / = (ii, . . . ,ik) the weight of I is 
|/| = ii + ■ ■ ■ + ik, and k = £{I) is the length of I. For / = {ii, . . . ,ik), the reversed 
composition {i^, ■ ■ ■ will be denoted J: of course reversal preserves weight, length and 
refinement (i.e., I J implies I J). We write J U J for the juxtaposition of / and J, 
so lUJ = Jul. 

Compositions of weight n are in 1-to-l correspondence with subsets of {1, 2, . . . , — 1} 
via partial sums 

{ii,i2,...,ik)^{ii,ii + i2,---,ii^ h ik-i}, (9) 

and this correspondence is an isomorphism of posets, i.e., J ^ J if and only if the subset 
corresponding to / contains that corresponding to J. Complementation in the power set 
thus gives rise to an involution / —>■ J*; e.g., (1,1,2)* = (3,1). Evidently |/*| = |/| and 
£{I) +£{!*) = |/| + 1. Also, J ^ J if and only if /* >z J*. The complementation operation 
commutes with reversal, so the notation I* is unambiguous. 

Because of the correspondence ([9]) between the poset of compositions of n and the 
poset of subsets of {1, 2, . . . , tt, — 1}, it follows that the Mobius function for compositions 
of n is given by 
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Thus, e.g., the Mobius inversion formula apphed to equation ([8]) is 

Mi = (10) 

Since QSym is a commutative Hopf algebra, its antipode S is an automorphism of 
QSym and S"^ = id. Now S can be given by the following explicit formulas: for proof see 
[5] or [15]. 

Theorem 3.1. The antipode S o/QSym is given by 

1. SiMj) = E/,u/,u...u/,=/(-l)'M,,M,, ■ ..Mj- 

2. S{Mj) = i-iy^'^Ej. 

Part (2) of this result implies a number of facts about the Ej. First, the Ej have almost 
the same multiplication rules as the Mj: if T is the automorphism of QSym sending Mj 
to Mj, then ST takes any identity among the Mj to an identity among the Ej that differs 
only in signs. For example, since 

M(2)M(3) = M(2,3) + M(3,2) + M(5) 

we have 

E{i)E(^i) = £'(2,3) + -£(3,2) — -£(5)- 

Second, the coproduct formula ([6]), which in terms of the Mj says 

A(M/) = ^ Mj® Mk, 

JUK=I 

can also be written 

A{Ej)= J2 Ej®Ek, 

JUK=I 

because of the standard Hopf algebra relation AS = {S ® S)A"p. Finally, if we apply S 
to both sides of Theorem 2.3 we obtain, for any composition / = (zi, . . . ,ik), 

Y,E..i= c{'B)M^,,^M^,,yM^,^^, (11) 

(TgSj. partitions S = {Bi, . . . , B;} of {1, . . . , fe} 

where as above c(!B) = (card-Bi — 1)! ■ ■ ■ (cardi?; — 1)! and bs = Tlji^Bs 
Now define an automorphism ip of 7j{x,y) by 

ll){x)=X + y, 'ijj{y) = -y. (12) 

Evidently tp"^ = id, and iplSj^) = S)^. Thus ip defines an additive involution of Sj^ = QSym 
(which is not, however, a homomorphism for the *-product). We can describe the action 
of ip on the integral bases for QSym as follows. 
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Theorem 3.2. For any composition I, 

1. ij{Mj) = 

2. ^{Ei) = -El*. 

Proof. Suppose w = is the word in x and y corresponding to a composition /. Then 
evidently substituting y in place of any particular factor x in w corresponds to splitting a 
part of I. With this observation, part (1) is clear (there is also one factor of —1 for each 
occurrence of y in w). 

Now we prove part (2). We have 

from part (1). From Example 1 of [15], S{Ff) = (—ly^^Fj*, where S is the antipode of 
QSym. Thus 

Now by Mobius inversion of equation ([7]) 
and so 

S^(Ej) = -(-l)^(^"*)Mj. 

Apply 5* be both sides to get 

□ 



4 Finite Multiple Sums and Mod p Results 

In this section we consider the finite sums 
and 



11 • • • n 

ri>ni>n2>--->nfc>l 1 k 



ri>ni>n2>--->nfc>l ^1 ' ' ' 



where the notation is patterned after that of [13]; the multiple zeta values are given by 

C(«i, ...,ik)= lim 

n— ►oo 

when the limit exists (i.e., when ii > 1). 

Let pn = eT(j)n, where 0„ is the map defined in §2, T is the automorphism of QSym 
sending Mj to Mj, and e is the function that sends tj to i. Then p„ : {S}^,*) — > R is 
a homomorphism sending Zi-^Zi,^ ' ' ' We can combine the homo- 

morphisms pn into a homomorphism p that sends w & Sj^ to the real- valued sequence 
n — >■ Pn{w). We shall write A/ for the real- valued sequence n — >■ ^/(n) (and similarly for 
S"/), so p0^^ sends M/ to A/ and Ei to S"/. For example, we can apply p0^^ to equation 
f fTOj) above to get 

A, = ^(-l)^(^)-^(^)^, (13) 

Applying the homomorphism Pn4'~^ to Theorem 2.3 and equation ( ITTi) . we get formulas 
for symmetric sums of A/(n) and Si{n) in terms of length one sums S'(m)(n) (cf. Theorems 
2.1 and 2.2 of [131). 



Theorem 4.1. For any composition I = {ii, ■ ■ ■ ,ik) , 

S^.i{n) = c(S)%)(n) ■ • • W 

o-eS;; partitions 2 = {Bi, . . . , _B;} o/ {1, . . . , fc} 

E ^-^H = H (-l)^"'c(B)%)(n) ■ ■ ■%)(n), 

(tSSj. partitions 2 = {Bi, . . . , -B;} o/ {1, . . . , fc} 



where c(!B) = (card-Bi — 1)! ■ ■ ■ (card 5; — 1)! and bg = Ylij^^Bs "3 



We consider two operators on the space of real- valued sequences. First, there is 
the partial-sum operator S, given by 



T,a{n) = a(z) 

i=0 

for a G R'^. Second, there is the operator V given by 

i=0 



It is easy to show that S and V generate a dihedral group within the automorphisms of 
RN, i.e., = id and SV = VS'^ It follows that (SV)^ = id. We have the following 
result on multiple sums. 

Theorem 4.2. For any composition I, SVS/ = —Sj*. 
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Proof. We proceed by induction on |/|. The weight one case is SV5'(i) = VS ^^(i) 
-5(1), i.e.. 



k \k J ^ k 

k=l ^ ' k=\ 



a classical (but often rediscovered) formula which actually goes back to Euler [6]. For 
i = (^i;^2; • • • j'^fc); it IS stralghtforward to show that VSj{n) = -V/(ra), where / G 
is given by 

S{i2,...,tk)in), ifii = l; 



fin) 



S ^'S'(ii_i,i2,...,ij^)(n), otherwise. 



Now suppose the theorem has been proved for all / of weight less than n, and let 
/ = (ii, . . . ,ik) have weight n. There are two cases: ii = 1, and ii > 1. In the first case, 
let (^2, • • • , ik)* = J = (ji, . . . , jr)- By the assertion of the preceding paragraph and the 
induction hypothesis, 

SV5,H = S(iv5j*(n)) = S(^S-15jH) = -%+i,,,,...,,„)(n). 

But evidently /* = (ji + 1, j2, • • • , jr), so the theorem holds in this case. 
If ii > 1, we instead write {ii — 1, ^2, • • • , h)* = J = {ji, ■ ■ ■ ,jr)- Then 

But in this case /* = (1, ji, . . . ,jr), so the theorem holds in this well. □ 

The proof of the preceding result is essentially a formalization of the procedure in 
App. B of [26]. Recalling the automorphism of S)^ defined by equation f|T2l) . we note 
that Theorem 4.2, together with part (2) of Theorem 3.2, says that the diagram 



QSym — QSym 



p 

rn rn 



(14) 



commutes. 

We now turn to mod p results about Sj{p — 1) and Af{p — 1), where p is a prime. Since 
the sums Aj{p — 1) and Si{p — 1) contain no factors of p in the denominators, they can be 
regarded as elements of the field Z/pZ. The following result about length one harmonic 
sums is well known (cf. [TTj, pp. 86-88). 

Theorem 4.3. S(^k)ip — 1) = mod p for all prime p > k + 1. 
Proof. Since the group of units in Z/pZ has order p — 1, we have 

xP-^ _ 1 = (a; _ l)(a; _ 2) ■ ■ ■ (x - p + 1) 
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in Z/pZ. Thus, in the expansion 

p 

x{x-l)---{x-p+l) = ^(-1)P~V41, 2, . . . ,p - l)x' 

1=1 



where Cj is the ith elementary symmetric function, we must have ep_j(l, . . . ,J9 — 1) = 
mod p except for eo(l, . . . ,p — 1) = 1 and ep_i(l, . . . ,p — 1) = —1 mod p. Thus, the 
Newton formulas 

T 

• • • ,P - l)er-i(l, . . . ,p - 1) = re^(l, . . . ,p - 1), 

i=l 

where pj is the zth power-sum symmetric function, reduce mod p to 

. . . ,p — 1) = — rer(l, . . . ,p — 1) = mod p 

when p > r + 1. But 

Pr(l, . . . ,p - 1) = 5'(r)(p - 1) modp 
since inversion permutes the units of Z/pZ, so the result follows. □ 

Combining the preceding result with Theorem 4.1 gives the following. 

Theorem 4.4. For any composition I = {ii, . . . ,ik) and prime p > |/| + 1, 

5^ A^.,(p-1)= 5^5^.,(p-l) = modp. 

In particular it follows that, for / = (r, r, . . . , r) {k repetitions), we have 

Ai{p-1) = Si{p-1) = modp (15) 

for prime p > rfc + 1 (cf. Theorem 2.14 of [2^). There is the following result relating 
sums associated to / and / (cf. Lemma 3.3 of [29]). 

Theorem 4.5. For any composition I , Aj(p—1) = (— l)l'^ly4/(p — 1) mod p, and similarly 
Sj{p — 1) = (— l)l'^l5'/(p — 1) mod p. 

Proof. Let I = {ii, . . . ,ik). Working mod p, we have 

0<bi<-<bk<p 1 k 

and similarly for Sj. □ 

11 



An immediate consequence is that Si{p — 1) = Aj{p — 1) = mod p if I = I and 
|/| is odd. Another consequence is that S(^ij){p — 1) = y4(jj)(p — 1) = mod p when 
p > i + i + 1 and i + j is even. This is because 

S(ij){p - 1) + -1)^0 mod p 

for p > 2 + j + 1 by Theorem 4.1, while S(^ij){p — 1) = SQ^i){p — 1) mod p when i + j is 
even by Theorem 4.5. 

We have the following result relating Sj and Si* . 

Theorem 4.6. Si{p — 1) = —Si*{p — 1) mod p for all primes p. 

Proof. Let / be a sequence. From the definition of V 

^v/H = X:(^|i) (-1)7(0, 

so taking n = p — 1 gives 

SV/(p - 1) = i-iy-'fip -l) = f{p - 1) mod p. 

Now take f = Sj and apply Theorem 4.2. □ 

Define, for each prime p, a map Xp : i^^ — > Z/pZ by Xp(^) = Pp-i('W^)- Then we can 
use the commutative diagram (JT^ to restate Theorem 4.6 in the following "algebraic" 
form, which corresponds to the duality theorem for multiple zeta values as formulated in 
Corollary 6.2 of [14j. 

Theorem 4.7. ^45 elements ofZjpZ, Xp{w) = Xp(V'(w)) for words w of . 
For example, since %p{x'^y^) = —x^y^ — xy'^ — yxy^ — y^, we have 

^(3,i,i)(p-l) = -^(3,i,i)(p-l)-A(2,i,i,i)(j9-l)-A(i,2,i,i)(p-l)-^(i,i,i,i,i)(p-l) mod p. 

(In fact, for p > 6 it follows from Theorem 7.1 below that both sides are congruent mod 
p to ^Bp^5, where Bi is the ith Bernoulli number.) 

5 Harmonic Sums of Height One 

Following the terminology of we say a word of S)^ of the form x^^^y^ has height one. 
We shall call the corresponding multiple harmonic sums sums of height one. 

Harmonic sums of this form have a number of special properties, which we discuss in this 
section. 

In this section only we use superscripts for repetition in compositions, so (n, l'^) means 
the composition of weight n + k with k repetitions of 1. For compositions of this form, 
we have the following result relating the two kinds of harmonic sums mod p. 
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Theorem 5.1. For any prime p > k, 

S(h,i^-i){p -l) = {-!)'' A(^h,ik-i){p - 1) mod p. 
Proof. Equate parts (1) and (2) of Theorem 3.1 and apply the homomorphism pp_i to get 

{-iY^'^Sj{p-i)= Yl (16) 

hU-Uli=I 

for any p. Now set / = {l''~^,h) in equation (fT6l) and reduce both sides mod p. The 
hypothesis insures that all terms on the right-hand side are zero mod p except the one 
with / = 1, and so we have 

(-l)^''^%^ifc-i)(p - 1) = A(^-^k-i^f,){p - 1) mod p, 

and the conclusion follows by Theorem 4.5. □ 

Combining the preceding result with the duality theorem gives the following congru- 
ence. It can be compared with Theorem 4.4 of [13], which asserts that ({h + 1, l'^"^) = 
({k + 1, l'*~^) for positive integers h and k. 

Theorem 5.2. If p is a prime with p > max{A;, h}, then 

\h,i^-^){p - 1) = A(^k,ih-i){p - 1) mod p. 
Proof. First note that {h, l^^'^y = (1^^^, k). So, combining Theorems 4.6 and 4.5, 

1) = 1) = 1) mod p. (17) 

Now by the preceding result 

Sih,i^-i){p - 1) = Af^h^ik-ijip - 1) mod p, 
- 1) = (-l)''^(fc,i'^-i)(p - 1) mod p, 

so the conclusion follows from congruence (ITTl) . □ 



The height one harmonic sums can be written in terms of Stirling numbers. The 
[unsigned) Stirling numbers of the first kind (Stirling cycle numbers) are given by 



number of permutations of {1, 2, . . . , n} with k disjoint cycles, 
and the Stirling number of the second kind (Stirling subset numbers) are given by 

Tl I 

^ > = number of partitions of {1, 2, . . . , n} with k blocks. 
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For basic properties of Stirling numbers see, e.g., [10]. Note that 



n + 1 
k+1 



e„_fc(l,2, . . . ,n). 



since each side is the coefficient of x'' in the expansion of x{x + 1) ■ ■ ■ (x + n). Hence 



n+1 
k + 1 



n>ji>--->j„_fc>l 



Jl • • -Jn-k 



n>ii>->ife>l 



n!A(ifc)(n). (18) 



An immediate consequence is the following result, which expresses height one sums in 
terms of Stirling numbers of the first kind. 

Theorem 5.3. For positive integers h, k, and n, 

m 



j=k 



Proof. We have 



E 



i^io ■ ■ - Ik — ' 1 

n>ii> ->ik>l ^ j=k 



^=i:iv-.,o-i)=E-i# 



j=k 



(where we used equation (|T8l) in the last step), from which the conclusion follows. □ 

In doing computations, the following mod p result is often useful. 
Theorem 5.4. For prime p and positive integers h,k < p, 

p—k 



1=1 K 3 ) 



p. 



Proof. We start by setting n = p—1 and reversing the order of summation in the preceding 
result: 

p—k 



By Fermat's theorem, [p — = {p — JY'^ = {—iY~-^ mod p. Also, since 



-1 



{p-l)\ 



p—1 

j-lj - (j-l)!(p-j)! 



mod p, 



it follows that 



{p-jV- 



{—iy{j — 1)! mod p 



14 



since (p — 1)! = — 1 mod p by Wilson's theorem. Hence 

p—k 

A(,,,.-.)(p-i) = ^(-iy(j-i)!(-jr 

i=i 

The conclusion then follows from congruence f|T9l) and a mod p relation between the two 
kinds of Stirling numbers: 

_jp-k\ 
= < > mod p 

IP — nj 

for prime p and l<k<n<p — 1, which follows by induction on n using the recurrence 
relations for the two kinds of Stirling numbers. □ 



P-J 
k 



mod p. 



(19) 



n 
k 



6 Results for Short Lengths 

In this section we prove some results about the sums Si{p — 1) for £{I) < 5, which will be 
useful for the calculations of the next section. (Similar results hold for the Ai{p — 1), but 
the Si{p — 1) are more convenient in view of Theorem 4.6.) Henceforth we will assume 
that all congruences are mod p. The following result expresses Sj{p — 1) in terms of a 
Bernoulli number when = 2 (cf. [29], Theorem 3.1). 

Theorem 6.1. For i,j positive and prime p > i + j + 1, 

Proof. We use the standard identity expressing sums of powers in terms of Bernoulli 
numbers (see, e.g., p!9]): 

a=l k=0 ^ ^ 

Using Fermat's theorem and equation fl20l) . we have 

p— 1 ^ a— 1 , p— 1 _ a— 1 

^,,)(p - 1) ^ A,,)(p -i) = E^E^-E^E 

a=l 6=1 a=l b=l 

p-1 1 T p-i-j / A 1 P-l-i / .\ P-1 



Now ^^ij a^~*~-'^^ = unless k = p — i — j . So the sum reduces to 



p-3\p-i-jj J\^J 
from which the conclusion follows. □ 
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Since odd Bernoulli numbers are zero, this result implies our earlier observation that 
S(^ij){p — 1) = when i + j is even. Note also that it implies S'(j^i)(p — 1) = for i 

even, so we can restate it as 

i.e., all double sums of odd weight n can be written in terms of 5'(„_i^i)(j9 — 1). In fact, 
if the weight is odd we can write all triple sums in terms of the same quantity (cf. [29j, 
Theorem 3.5). 

Theorem 6.2. Ifn = i + j + k is odd and p > n + 1, then 



%i,fc)(p-l) = ^ 



S{n-l,l){P - !)• 



Proof. From Theorem 3.1 it follows that 

Mm)= i-lf'^-'Si^im)Sj^{m)...Sj^im) (22) 

IlUl2U---Uli=I 

for any composition / and positive integer m. In particular, for £(/) = 3 we have 

A^ij^k){m) = S'(i)(m)S'(j)(m)S'(fc)(m) - % i)(m)S'(fc)(m) - S^i){m)Si^k,j){m) + S(^k,j,i){m). 
But also, from equation (fT3|) . 

-4(ij,fe)(m) = S(_ij^k){m) - 5'(i+j-fc)(m) - 5'(ij+fc)(m) + S(^i+j+k){m). 
In the case m = p — 1, this gives the congruence 

S{k,j,i)iP - 1) = 'S'(ij- fe)(p - 1) - 5'(i+j-fc)(p - 1) - S(^ij+k)ip - 1), 
using Theorem 4.3, and if in addition i + j + k is odd, we have 

2^(i,j,fc)(p - 1) = S(^i+j^k){p - 1) + S(^ij+k){p - 1) 
from Theorem 4.5. Now use congruence (12T!) . □ 

For a composition / = {ii, i2, . . . , ik) of odd weight and length > 1 it is convenient 
to define 

fc-l ^ I ^1 X 

(23) 



j=l ^ 



so Theorem 6.2 says that 



5/(p-l) = ^C^(/)%-i,i)(p-l) 



when n = \I\ is odd and = 3. It is easy to show that the function C{I) mirrors 
the properties of Si{p — 1) with regard to reversal and duality, i.e., C(/) = — C(/) and 
C(/*) = -C{I). 

In the even-weight case, we can express all length-4 sums in terms of triple sums and 
products of lower-weight sums as follows. 
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Theorem 6.3. Ifn = i + j + k + l is even and p > n + 1, then 

2S(^ij^k,i){p - 1) = S(i+j^k,i)ip - 1) + S(^ij+k,i){p - 1) + 

S(i,j,k+i){p - 1) + S(^ij){p - - 1). 

Proof. From equation (I22p with £{I) = 4, 

A(ij^k,i){rn) = -S(i^k,j,i){m) + S(i){m)S(i.j^k){m) + S(j^i){ni)S(i^k){'m) + S(^k,j,i){m)S(k){m) 

- S'(i)(m)S'(j-)(m)S'(;,fc)(m) - S'(i)(m)S'(fej)(m)S'(,)(m) - % i)(m)S'(fc)(m)S'(i)(m) 

+ S^i){m)S^j){m)S(^k){m)S^i){m). 

Set m = p — 1 and reduce mod p to get 

^(i,i,fc,/)(P - 1) = -S{l,k,j,i)iP - 1) + - l)'S'(Z,fc)(p - 1). 

But equation f|T3l) gives, after reduction mod p, 

^(ij,fc,o(P - 1) = S(^,j^k,i){p - 1) - - 1) - S(^ij+k,i){p - 1) - S(^ij^k+i){p - 1), 

and the conclusion fohows by the use of Theorem 4.5. □ 
If we speciahze the preceding resuh to the case k = i,l = j, we have 

2%,j,ij)(P - 1) = - 1) + - 1) + - 1) + - 1)^- 

But the sum of the first three terms on the right-hand side is zero mod p by Theorem 4.4, 
so we have 

2%,,,,)(p-l) = %,)(p-l)2. (24) 

We note that we also have A(^ij^ij){p — 1) = S(ij^ij){p — 1), which follows from equation 
( IT3ll . Theorem 4.4, and the comments following Theorem 4.5. If i and j have the same 
parity, then i + j is even and the right-hand side of (12^ is zero mod p by Theorem 4.5, 
so in this case 

Sii,j,i,j)iP - 1) = ^{i,j,i,j)iP -1)^0. 

Cf. Theorem 3.18 of [29j. 

It is evident that the ideas of the two preceding results can be used to write Sj{p — 1), 
where |/| and are of the same parity, in terms of sums of shorter lengths, together 
with products of sums of lower weight. We state one more result of this type, since we 
need it for our computations in weight 9; we omit the proof, which is straightforward. 

Theorem 6.4. If I = {i, j, k, I, m) is a composition of odd weight n and p > n + 1, then 

2Sl{p — 1) = S(iJ^j^k,l,m){P — 1) + S(ij+k,l,m){P — 1) + S(ij^k+l,m){p " 1) + S(ij^k,l+m){p " 1) 

- ^C{I)S(^n-l,l){P - 1) + S(^ij){p - l)S(^k,l,m){p - 1) + S(ij^k){P - l)^{«,m)(p - 1), 

where C{I) is given by equation l[2^). 
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7 Calculations in Low Weights 



In this section, we find, for n < 9, sets of quantities tliat generate all multiple harmonic 
sums Si{p — 1) of weight n mod p when p > n + 1. In fact, all our generators will be 
sums of the form S(2h,i,...,i){p — 1) (= A{2h,i,...,i){p — 1) by Theorem 5.1), or products of 
such sums. To make the notation less cumbersome we will use an abbreviated form in 
referring to specific sums, e.g., we write 5*221 instead of S(^2,2,i){p — 1)- We will also assume 
that p > |/| + 1 in considering any sum Si{p — 1). 

The congruence (fT5|) implies Si{p — 1) = when / has length 1 or |/|, so the only 
sums Si{p — 1) with |/| < 3 that are (possibly) nonzero mod p are 5*21 = -Bp-s and 
5*12 = —Bp^^. Now p divides the numerator of -Bp-3 only for certain irregular primes p, 
the only examples with p < 1.2 x 10^ being p = 16, 843 and p = 2, 124, 679 [iHl H [3] . 

For / of weight 4 we have Sj{p — 1) = 0; this follows from congruence f|T5|) if = 1 
or 4, by the remarks following Theorem 4.5 if = 2, and by duality if = 3 (since 
in that case £{!*) = 2). 

In weight 5, it suffices by duality to consider Si{p — 1) with < 3. Theorems 6.1 
and 6.2 give the following result, which implies that all weight 5 sums vanish mod p if 
and only if p divides the numerator of -Bp-5; the only known p for which this happens is 
p = 37, and this is the only such p < 1.2 x 10^ p]. 

Theorem 7.1. All Sj{p — 1) with \I\ = 5 are multiples of Sn = Bp_^. In particular, 
S32 = —25*41, 5311 = — |54i, and S221 = |54i. 

For |/| = 6, it follows from the remarks after Theorem 4.5 and duality that 5/(p— 1) = 
unless £(/) = 3 or 4. By duality it suffices to consider length 3. Here we have the following 
result , which implies that all harmonic sums of weight 6 vanish if and only if p divides 

Bp-3- 

Theorem 7.2. All the Si{p — 1) with \I\ = 6 are multiples of 5*4ii = —^Bp_^. In 
particular, 5*i4i = — 25*4ii, 5*312 = — 5*4ii, 5*321 = — 25*4ii, and 5*231 = 35*4ii. 

Proof. From Theorem 4.4, the symmetric sum 

5411 + 5i4i + 5ii4 = 25411 + 5i4i 

is congruent to zero, from which the statement about 5*i4i follows. For the same reason 

53111 + 5i3ii + 5ii3i + 5iii3 = 253111 + 25i3ii 

is congruent to zero; applying duality, we get 5*4ii + 5*312 = 0. Theorem 4.4 also gives 

5*321 + 5*231 = —5*312 = 5411 (25) 

Now multiplying in QSym and reducing mod p, we have 

=5*i5*22i 

= 5l221 + 52121 + 252211 — 5*321 — 5*231 — 5*222 
= 5*2121 + 25*2211 — 5*321 — 5*231. 
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Apply duality (and Theorem 4.5) to obtain from this 

3^321 + 2^231 = 0, (26) 



and solve congruences fl25|) and fl26l) to obtain 6*321 and 6*231 in terms of 6*411. 

Finally, 6*2121 = ^^21 by congruence fl24l) . while from duality we have 62121 = —5*231 



Since we've already shown 6*231 = 36411, the congruence S21 = —66411 follows. □ 



5*52 




— 35*61 


643 




5661 


5511 




— 561 


5*421 




35*61 


6412 




661 


5241 




2661 


5*331 




— 25*61 


5322 




-4661 



In weight 7, we need only consider sums of length 4 or less. By Theorems 6.1 and 6.2 
we can express all sums of length 2 and 3 in terms of 5*6i = Bp_j, as follows. 

(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 

For the length 4 sums we have the following result. 

Theorem 7.3. The length 4 harmonic sums of weight 7 can be written in terms of Sqi = 
Bp^-j and 64111. In particular 63121 = 62221 = 64111, and 

5*1411 = 3661 — 364111 
5*1321 = 95*61 — 35*4111 
5*1312 = —95*61 + 55*4111 
63211 = —6661 + 64111. 

Proof. Duality gives 6*1411 = 5*3ii2, 6*3121 = 6*2311, and 6*2212 = 5*i32i. Thus, it suffices 
to consider the seven length-four sums 6*4111, 5*i4ii, 6*3121, 6*2212, 5*i3i2, 6*2221, and 6*3211. 
We can obtain relations among these sums and 6*51 by multiplying lower-weight sums in 
QSym, reducing mod p, and then applying congruences (27-34) above to get all terms of 
length three or less in terms of 6*51. For example, multiplying 5*i by 6*411 gives 

35*61 = 5*1411 + 35*4111. (35) 

Similarly, consideration of the products 626311, 626221, 6*21631, and 6*11632 gives respec- 
tively 

I- 5*3211 + 5*1411 (36) 

^ 5*2212 (37) 

+ 35*3121 + 263211 (38) 

5*2212 + 5*1312 + 5*3211. (39) 



—35*61 


^ 25*3121 - 


9^61 


= 35*2221 - 


—35*61 


= — 5i3i2 


—35*61 


= 5*1411 + 
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Finally, multiplying 5*1 by S^m and using duality gives 

— 3Sqi = S'4111 + S'i4ii + S'3211 + 5'3i2i (40) 

The six congruences (35-40) can then be solved for Sqi and S'4111 to obtain the conclusion. 

□ 

Remark. While it does not appear that S'4111 is a multiple of Sei, we have not been able 
to prove this. There is no known prime p for which -Bp-y = 0; in particular, Bp^j ^ for 
allp < 1.2 X 10^ [3]. 

In weight 8, it suffices by duality to consider / with i{I) < 4. By Theorem 6.3, we 
can write all length 4 sums in terms of triple sums and S4iS'2i, and of course all double 
sums are zero. So it is enough to consider length 3, where we have the following result. 

Theorem 7.4. All length 3, weight 8 sums can be written in terms of Squ and 84^1821 = 
Bp^^Bps. In particular, 

S'521 = S'eii + S'4iS'2i (41) 

S'512 = -{—TSqu — S'4iS'2i) (42) 

^431 = ^(-25^611 -9^41^21) (43) 
S413 = ^(5S611 + S4lS2l) (44) 

S'332 = S'422 = lOS'eii + 2S4iS'2i. (45) 

Proof. Consideration of the products S'iS'211111 = and S'iS'121111 = 0, using duality, gives 
respectively 

Sen + S521 + S'431 + S341 + S251 + S161 = 0, (46) 
S'521 + S'512 + S'422 + S'332 + S'242 + S'i52 = 0. (47) 

Now from symmetry (Theorem 4.4), S'lgi = — 2S'6ii, S'251 = —S'521 — S'512, and S'341 = 
—S'431 — S'413. Using these facts, congruences (46-47) reduce to 

S'413 = Sen — S512, (48) 

^332 = S422. (49) 

Now consider the product S'iS'31111 = 0. Using duality, this gives 

5S611 + S512 = S4211 + S3311 + S2311. 

Using Theorem 6.3, the right-hand side can be expressed in terms of triple sums and 
S'4iS'2i, and then applying the preceding relations we obtain 

S'521 = — 6S'6ii — 2S'5i2. (50) 
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Applying similar techniques to the product <S'25'2iii and using congruences (48-50) gives 
congruences (41), (42), and (44). 

Finally, using duality on the product 5'25'2ii2 = gives 

= 4jS'i42i + 2jS'4ii2 + 2jS'2312 + 'S'laai- 

Now rewrite this using Theorem 6.3 and apply the previously established relations to 
obtain (43); congruence (45) then follows. □ 

Remark. It cannot be the case that Sen — US41S21, where w is a unit of Z/pZ for all 
p > 9, since when p = 37 we have Squ = 7 and 5'4i = 0. 

For sums of weight 9, duality allows us to restrict attention to 5'/(p — 1) with £{I) < 5. 
Using Theorem 6.4 (and Theorem 7.2), we can rewrite the length 5 sums in terms of length 
4 sums, Ssi, and the product 5'4iiS'2i. So it is enough to consider the sums Si{p — 1) with 
= 4, for which wc have the following result. 

Theorem 7.5. All length 4, weight 9 harmonic sums can he written in terms of the three 
quantities Sqiu, S^i = -Bp-g, and 5'2i5'4ii = —^B^_^. In particular, 

19 310 1 

'S'leii = ~35'6iii + "^'S'si 'S'5211 = <S'6iii — 27"*^^^ 2>S'2i<S'4ii 

'S'5112 = — 45'6iii + ^^'^si — g'S'215'411 >S'5i2i = 5'6iii — 27*^*^1 ^ 25'2i5'4ii 

89 1 544 1 

'S'2511 = 25'6iii — ^'S'si — -5'2iS'4ii 'S'2151 = — 75'6iii + ^^"^si ~ 2*S'2iS'4ii 

137 13 

'S'1521 = — 35'6111 H g~'^81 ~ 'S'21'S'411 'S'4311 = "S'eiii + "g"'^81 

56 133 

'S'4131 = "65*6111 + "^'S'si — >S'2i5'4ii 54113 = 7^6111 g"*^^-^ 'S'215411 

20 182 

53411 = —256111 — g"*^^^ 'S'3141 = 856111 g"*^^-^ 'S'2i54ii 

5i43i = 456111 3~'^S1 52l54ii 54221 = 56111 + ^'S'si — -52l54ii 

617 2 370 1 

54212 = —456111 + "^^9,1 ~ -52i54ii 54122 = 356111 — 27"*^^^ 2'S'2i54ii 

124 4 110 

52421 = 256111 2^*^^^ g'S'2l54ii 52412 = "256111 H g^'^Sl 

_ 122 2 _ 290 2 

52241 = "27"*^^^ ~ 2'S'2i54ii 53321 = —256111 + ^^'S'si ~ ^5215411 

289 1 292 4 

53312 = 256111 — 27"*^^^ 2'S'2i54ii 53231 = ""27"*^^^ ^ 2>S'2i54ii 

_ 1 _ 64 

52331 = g'S'si ~ 52l54ii 53213 = 456111 3"*^^^ 

53132 = "27'^^''^ ~ 2'S'2l54ii 52322 = — -g-581 

53222 = —6581. 
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Proof. There are 56 length-4 sums, but by reversal (Theorem 4.5) we can reduce to 28 
of them. Next we go from 28 lcngth-4 sums to 14 length-4 sums, together with Sgi, by 
considering various products in QSym. For example, from the products 5'i 5*422, 'S'25'421, 
and 8284^12 respectively, we obtain the relations 

23 

-^'S'si = —82241 + 'S'4122 + 54212 + ^4221 

73 

— 581 = 52421 + 254221 + 54212 

23 

-g-581 = 52412 + 54212 + 254122 

which we can use to write 5224i, 5242i, and 52412 in terms of 54221, 54212, 54122, and 58i. 
Similarly, we eliminate 5i6ii, 525ii, 52i5i, 5i52i, 53411, 53141, 5i43i, 52331, 53213, 53132, 
and 52322- Now we consider the eight products 5i53iiiii, 5i522iiii, 5i5i3iiii, 5i5ii3iii, 
5i 5212111, 5i 5211211, 5i 5211121, and 5i5i22iii- Using duality and rewriting in terms of our 
14 length-4 sums, we have the relations 

28 

-^581 = 256111 + 55121 + 55112 + 54131 -|- 54113 

-^581 = 55211 + 54221 + 54212 ~ 53321 
29 

-^581 = —55112 — 254122 + 53312 
6 1 

— -^581 = —54113 + 54122 + 253312 

4 _ 

^581 = 54311 — 54221 + 53321 + 53312 

11581 = —54113 + 53141 + 54311 -|- 53321 + 53231 + 54122 + 54212 + 54221 

—13581 = 55211 + 55112 + 55121 + 54122 — 54221 — 53231 — 53321 — 53312 

— 6581 = 53222- 

We can use these congruences to solve for all the length-4 sums in terms of the seven 
quantities 58i, 56iii, 55211, 55112, 55121, 54311, and 54221- 

Prom duality, 541211 = 533111. Expanding out both sides using Theorem 6.4, we obtain 
the congruence 

55211 + 54311 -|- 54131 -|- 54122 = 56111 + 53411 -|- 53321 + 53312, 

which (after rewriting in terms of our seven chosen quantities), allows us to solve for 54221 
in terms of the remaining six quantities. Now apply Theorem 6.4 similarly to the duality 
relations 541112 = 5i5iii, 5i24ii = 531122, and 52i4ii = 531131 to obtain respectively 

—54311 -|- 255121 + 55211 — 2561111 + 21581 = 52i54ii (51) 

91 

354311 — 55121 + 455211 — 6561111 + "^'S'si = 52i54ii (52) 

89 

—354311 -|- 255121 + 55211 + "2"'^8i = 52i54ii. (53) 
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Add 3 times congruence (51) to congruence (52), and compare the resulting expression 
for 5*21 5*411 with that obtained by adding congruence (53) to (52); the result is 

80 

>S'5l21 = >S'5211 + —Ssi. (54) 

y 

Next substitute congruence fl5^ into the result of subtracting congruence (52) from (51) 
to get 

13 

5'43ii = Squi + "^'S'si- (55) 

Using congruences and ([S5D, any of the congruences (51-53) can be used to write 
5*5211 in terms of 5*6111 , 5*81, and 5*215*411. 

Finally, consideration of the product 5*ii 5*2221 gives the relation 



35*3231 + 35*3321 + 35*4221 + 25*3222 + 25*2331 + 5*2241 + 5*2421 

= 3531221 + 3532121 + 25*32211 + 2522131 + 5*21231 + 5*23121, 

which can be used to write 55112 in terms of 56iii, 58i, and 52i54ii. □ 

Remarks. 1. The only primes p < 1.2 x 10^ with Bp_g = are p = 67 and p = 877 
[3]. Here is a table of the values of our generators at these primes, and at p = 16,843, 
computed using Mathematica and Theorem 5.4. 



p 


67 


877 


16,843 


5*6111 


7 


253 


16,690 


581 








14,820 


52l54ii 


4 


354 






2. J. Zhao conjectures that 



c _ 1 c c I ^^^^ c 

•Jeill = — -'-'21'-'411 + „.„ ^81 

9 648 



(see [2^ Prop. 3.12]). This conjecture is consistent with the data in the preceding remark. 
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